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Abstract. Reconstructions of Quaternary climate are of- 
ten based on the isotopic content of paleo-precipitation pre- 
served in proxy records. While many paleo-precipitation 
isotope records are available, few studies have synthe- 
sized these dispersed records to explore spatial patterns 
of late-glacial precipitation <$ 18 0. Here we present a syn- 
thesis of 86 globally distributed groundwater (n = 59), 
cave calcite (n = 15) and ice core (n = 12) isotope 
records spanning the late-glacial (defined as ~ 50 000 
to ~ 20 000 years ago) to the late-Holocene (within the 
past ~5000 years). We show that precipitation <$ 18 0 
changes from the late-glacial to the late-Holocene range 
from — 7.1 %c (5 18 Oi a t e -Holocene > <5 18 0i a t e -glacial) lo Tl.7 %c 
(5 18 Oi at e-giacial>5 18 Oi ate -Holocene), with the majority (77%) 
of records having lower late-glacial <5 18 0 than late-Holocene 
<5 18 O values. High-magnitude, negative precipitation <5 18 0 
shifts are common at high latitudes, high altitudes and con- 
tinental interiors (5 18 Oi ate -Holocene >^ 18 Oi a te-glacial by more 
than 3 %c). Conversely, low-magnitude, positive precipita- 
tion <$ 18 0 shifts are concentrated along tropical and sub- 
tropical coasts (5 1 8 Oi at e-gl a cial > <5 1 8 Oi ate -Holocene by less than 
2 %c). Broad, global patterns of late-glacial to late-Holocene 
precipitation <$ 18 0 shifts suggest that stronger- than-modem 
isotopic distillation of air masses prevailed during the late- 


glacial, likely impacted by larger global temperature differ- 
ences between the tropics and the poles. Further, to test how 
well general circulation models reproduce global precipita- 
tion 5 18 0 shifts, we compiled simulated precipitation <$ 18 0 
shifts from five isotope-enabled general circulation mod- 
els simulated under recent and last glacial maximum cli- 
mate states. Climate simulations generally show better inter- 
model and model-measurement agreement in temperate re- 
gions than in the tropics, highlighting a need for further re- 
search to better understand how inter-model spread in con- 
vective rainout, seawater <5 18 0 and glacial topography pa- 
rameterizations impact simulated precipitation <$ 18 0. Future 
research on paleo-precipitation <$ 18 0 records can use the 
global maps of measured and simulated late-glacial precipi- 
tation isotope compositions to target and prioritize field sites. 


1 Introduction 

Isotopic compositions of late-glacial precipitation can be pre- 
served in groundwaters, cave calcite, glacial ice, ground ice 
and lake sediments. These records have been used to better 
understand past climate changes for more than a half cen- 
tury (e.g., Munnich, 1957; Thatcher et al., 1961; Munnich et 
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al., 1967; Pearson and White, 1967; Tamers, 1967; Gat et al., 
1969). Each type of isotopic proxy record is distinguished by 
its temporal resolution, preservation of one or both 18 O / 16 O 
and 2 H / 1 H ratios, and frequency on land surface. For exam- 
ple, groundwater records contain both 18 0 / 16 0 and 2 H / 
ratios with widespread global occurrence, but have a coarser 
temporal resolution than other paleoclimate proxies (Rozan- 
ski, 1985; Edmunds and Milne, 2001; Edmunds, 2009; Cor- 
cho Alvarado et al., 2011; Jirakova et al., 2011). Speleothem 
records, by contrast, have high temporal resolution but usu- 
ally only report calcite 18 O / 16 O ratios (without fluid inclu- 
sion 2 H / 1 H data) and are less common than groundwater 
records (e.g., Harmon et al., 1978, 1979). Late-glacial ice 
core and ground ice records have high temporal resolution, 
can be analysed for 18 0 / 16 0 and 2 H / l H ratios, but are rare 
on non-polar lands (Dansgaard et al., 1982; Thompson et al., 
1989, 1995, 1997, 1998). Lake sediment records can have a 
high temporal resolution, can preserve 18 0 / 16 0 and 2 H / l H 
ratios and are available for a multitude of globally distributed 
locations (e.g., Edwards and McAndrews, 1989; Eawag et al., 
1992; Menking et al., 1997; Wolfe et al., 2000; Anderson et 
al., 2001; Beuning et al., 2002; Sachse et al., 2004; Morley 
et al., 2005; Tierney et al., 2008). However, some lake water 
proxy isotope records may be impacted by paleo-lake evap- 
orative isotope effects that obscure the primary meteoric wa- 
ter signal and mask paleo-precipitation isotope compositions 
(e.g., lake sediment calcite, diatom silica; Leng and Marshall, 
2004). 

This study examines speleothem, ice core and ground- 
water isotope records, focusing primarily on the groundwa- 
ter isotope records due to their relative density in the pub- 
lished literature in comparison to the more limited number 
of published speleothem and ice core records (compilations 
by Pedro et al., 2011; Stenni et al., 2011; Clark et al., 2012; 
Shah et al., 2013; Caley et al., 2014a). There exist roughly 
twice as many groundwater reconstructions of late-glacial to 
late-Holocene precipitation 5 18 O shifts ( n = 59) as the com- 
bined total of speleothem and ice core records (n = 27 ; where 

^ ^ O = (^O/ 1 ^ O sam pi e ) / (^O/ 1 ^ O standard mean ocean water" 

1) x 1000). A recent global synthesis of paired precipitation- 
groundwater isotopic data demonstrated that modern annual 
precipitation and modern groundwater isotope compositions 
follow systematic relationships with some bias toward win- 
ter and wet-season precipitation (Jasechko et al., 2014). Sys- 
tematic rainfall-recharge relationships shown by Jasechko et 
al. (2014) support our primary assumption in this study that 
groundwater isotope compositions closely reflect meteoric 
water. Because groundwater records can only identify cli- 
mate change occurring over thousands of years due to hydro- 
dynamic dispersion during multi-millennial residence times 
(e.g., Davison and Airey, 1982; Stute and Deak, 1989), we 
limit the focus of this study to meteoric water isotope compo- 
sition changes from the latter half of the last glacial time pe- 
riod to the late-Holocene. The latter half of the last glacial pe- 
riod is defined as ~ 20 000 to ~ 50 000 years before present, 
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using the end of the last glacial maximum as the more recent 
age limit (~ 20000 years before present; Clark et al., 2009) 
and the maximum age of groundwater that can be identified 
by 14 C dating as an approximate upper age limit (i.e., ground- 
water ages more recent than ~ 50 000 years old). 

For brevity, we refer herein to the time period repre- 
senting the latter half of the last glacial period (~ 20 000 
to ~ 50 000 years before present) as the late-glacial (e.g., 
5 1 8 Oiate-giaciai ) • We adopt a definition of the late-Holocene 
as occurring within the last 5000 years following Thomp- 
son et al. (2006). Other work proposes the late-Holocene be 
defined as within the last 4200 years (Walker et al., 2012), 
which is consistent with the 5000 years before present defini- 
tion (Thompson et al., 2006) within the practical uncertainty 
of 14 C-based groundwater ages (=b ~ 10 3 years). Further, al- 
though precipitation isotope compositions have varied over 
the late-Holocene, groundwater mixing integrates this vari- 
ability, prohibiting paleoclimate interpretation at finer tem- 
poral resolutions. 

Late-glacial to late-Holocene changes in precipitation iso- 
tope compositions provide important insights into conditions 
and processes of the past. Perhaps the two best-constrained 
global-in- scale differences between the late-glacial and the 
late-Holocene are changes to oceanic and atmospheric tem- 
peratures (MARGO Members, 2009; Shakun and Carlson, 
2010; Annan and Hargreaves, 2013), and changes to seawater 
<$ 18 0 (Emiliani, 1955; Dansgaard and Tauber, 1969; Schrag 
et al., 1996, 2002). Atmospheric temperatures have increased 
by a global average of ~ 4 °C since the last glacial maximum, 
with greatest warming at the poles and more modest warming 
at lower latitudes (Fig. 1; Shakun and Carlson, 2010; Annan 
and Hargreaves, 2013). Seawater <$ 18 0 during the last glacial 
maximum was 1.0 zb 0.1 %c higher than the modern ocean, 
as constrained by paleo-ocean water samples collected from 
pore waters trapped within sea floor sediments (Schrag et al., 
2002 ). 

Previous studies have proposed many different interpreta- 
tions of past changes to precipitation isotope compositions. 
Records of paleo-precipitation <S 18 0 have been used as a 
proxy for regional land surface and atmospheric tempera- 
ture (e.g., Rozanski, 1985; Nikolayev and Mikhalev, 1995; 
Johnsen et al., 2001; Grasby and Chen, 2005; Akouvi et 
al., 2008; Bakari et al., 2012); however, <$ 18 0-based pale- 
otemperatures can be complicated by past changes to a vari- 
ety of other processes controlling precipitation <$ 18 0, includ- 
ing moisture sources, upwind rainout, transport pathways, 
moisture recycling and in-cloud processes (Ciais and Jouzel, 
1994; Masson-Delmotte et al., 2005; Sjostrom and Welker, 
2009). Process-based explanations for observed meteoric wa- 
ter 8 18 O variations in proxy records include changes to hurri- 
cane intensity (e.g., Plummer, 1993), large-scale atmospheric 
circulation (e.g., Rozanski, 1985; Weyhenmeyer et al., 2000; 
McDermott et al., 2001; Pausata et al., 2009; Asmerom et al., 
2010; Oster et al., 2015), aridity (e.g., Wagner et al., 2010), 
monsoon strength (e.g., Denniston et al., 2000; Lachniet et 
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Last glacial maximum temperature minus modern (°C) 


Figure 1 . The change in surface air temperatures from the last glacial maximum to the preindustrial era (gridded data from Annan and 
Hargreaves, 2013). (a) Percentile ranges of temperature changes since the last glacial maximum for 10 degree latitudinal bands. Blue shading 
marks the 25th-75th percentile range; thin horizontal lines mark the 10th-90th percentile range. The grey band shows the globally averaged 
estimate of temperature change since the last glacial maximum of —4.0 =b 0.8 °C. (b) Gridded surface air temperature anomaly from the last 
glacial maximum to the preindustrial era (data from Annan and Hargreaves, 2013). 


al., 2004; Liu et al., 2007; Pausata et al., 2011a), local seawa- 
ter <$ 18 0 (Wood et al., 2003; Feng et al., 2014), precipitation 
seasonality (e.g., Fawcett et al., 1997; Werner et al., 2000; 
Cruz et al., 2005), moisture provenance (e.g., Sjostrom and 
Welker, 2009; Lewis et al., 2010), storm tracks, climate oscil- 
lation modes (e.g., North Atlantic oscillation), moisture recy- 
cling (e.g., Winnick et al., 2013, 2014; Liu et al., 2014a, b) 
and groundwater flow path architecture (Purdy et al., 1996; 
Stewart et al., 2004; Morrissey et al., 2010; Hagedorn, 2015). 
While unravelling these mechanisms and delineating the pri- 
mary and secondary processes can be rather challenging, the 
use of climate models in combination with robust and ex- 
tensive precipitation isotope data can resolve many of these 
complexities with meaningful interpretations and insight. 

The objective of this study is to analyse spatial patterns 
of measured late-glacial to late-Holocene precipitation <5 18 0 
changes from published groundwater, ground ice, glacial ice 
and cave calcite records, and to compare these measurements 
with output from five state-of-the-art isotope-enabled general 
circulation model simulations of last glacial maximum and 
pre-industrial or modern climate conditions. Synthesizing 
paleowater <$ 18 0 records provides an important constraint for 
isotope-enabled general circulation model simulations of at- 
mospheric and hydrologic conditions during glacial climate 
states (Jouzel et al., 2000). We combine a new global compi- 
lation of late-glacial groundwater and ground ice isotope data 
( n = 59) with existing compilations for speleothems (n = 15; 
Shah et al., 2013) and ice cores (n = 12; Pedro et al., 2011; 
Stenni et al., 2011; Clark et al., 2012; Caley et al., 2014a). 
This compilation of late-glacial groundwater isotope compo- 
sitions builds from earlier reviews of European and African 
paleowater isotope compositions (Rozanski, 1985; Edmunds 
and Milne, 2001; Darling, 2004; Edmunds, 2009; Negrel and 
Petelet-Giraud, 2011; Jirakova et al., 2011). 


2 Data set and methods 

In order to examine spatial patterns of change to mete- 
oric water <$ 18 0 values we compiled <$ 18 0, 8 2 H, 8 l3 C and 
14 C data from 1713 groundwater samples collected from 59 
aquifer systems reported in 76 publications (data and pri- 
mary references presented in the Supplement). <$ 13 C, 3 H and 
14 C data were used to estimate groundwater age (details 
within Supplement). Changes to precipitation <$ 18 0 values 
over time were determined by comparing groundwater iso- 
tope compositions of the late-Holocene (5 18 Oi ate -Hoiocene de- 
fined here as less than 5000 years before present; Thompson 
et al., 2006) and the latter half of the last glacial time period 
(5 1 8 Oi a te-giaciai ■ 20 000 to ~ 50 000 years before present). We 
acknowledge that these two relatively long time intervals - 
necessarily long in order to examine groundwater isotope 
records - integrate precipitation <$ 18 0 variability over the 
course of each time interval. The late-Holocene time interval 
integrates known precipitation <$ 18 0 variability (e.g., Aichner 
et al., 2015), and the late-glacial time interval likely incorpo- 
rates groundwater preceding the last glacial maximum, po- 
tentially during Marine Isotope Stage 3 or even older glacial 
time periods due to large uncertainties in 14 C-based ground- 
water ages (Supplement). 

Proxy-based meteoric water <$ 18 0 changes from the lat- 
ter half of the last glacial time period to the late-Holocene 
are described herein as measured A 18 Oi a t e -glaciah where 
measured A^Oiate-glacial = <^Oi a te-glacial — ^late-Holocene* 
A minimum groundwater age of 20 000 years before present 
was used to define the late-glacial to remain consistent with 
the timing of the last glacial maximum (~ 20 000 years be- 
fore present; Clark et al., 2009). Samples having a deu- 
terium excess of less than zero (deuterium excess = 8 2 H 
— 8 x <5 18 0; Dansgaard, 1964) and falling along region- 
ally characteristic evaporation <$ 2 H/<$ 18 0 slopes (Gibson et 
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al., 2008) were removed from the analysis to avoid includ- 
ing groundwater samples impacted by partial evaporation. 
Further, studies reporting saltwater intrusion were avoided 
on the basis of groundwater <5 18 0 and salinities showing 
evidence of seawater mixing (e.g., Schiavo et al., 2009; 
Yechieli et al., 2009; Hamouda et al., 2011; Han et al., 2011; 
Wang and Jiao, 2012; Currell et al., 2013). The 59 compiled 
groundwater measured A 18 Oi a t e -giacial values are unevenly 
distributed among western Europe (n = 10), eastern Europe 
and the Middle-East (n = 12), Africa (n = 11), southeastern 
Asia (n = 6), Australia, Oceania and the Malay Archipelago 
(n= 2), South America (n = 2), temperate and subtropical 
North America ( n = 8) and the High Arctic ( n =2). Half of 
the compiled groundwater records are located in the tropics 
or subtropics (that is, within 35° of the equator; n = 29) and 
half are located in the extra-tropics ( n = 30). 

Speleothem and ice core isotope proxy records were also 
compiled. Lacustrine sediment <$ 18 0 records are not consid- 
ered in this study because these records may preserve me- 
teoric waters impacted by evaporative isotope effects (Leng 
and Marshall, 2004). Speleothem and ice core measured 
A 1 8 Oiate-giaciai values were calculated by subtracting average 
<$ 18 0 values for each of the two time intervals defined for the 
groundwater records: the late-Holocene (< 5000 years before 
present) and latter half of the last glacial time period (20 000 
to 50 000 years before present). This step effectively lowered 
the temporal resolution of speleothem and ice core precipi- 
tation isotope records to be consistent with the temporal res- 
olution of the groundwater records. A correction factor was 
applied to speleothem <$ 18 0 values to account for different 
H 20 -CaC 03 isotopic fractionation factors during the late- 
glacial and the late-Holocene because of differing land sur- 
face temperatures during each time period (details presented 
within Supplement). 

Simulated A 18 0i a te-gladal values were compiled from 
five isotope-enabled general circulation models (simulated 

A Oi a te-glacial = ^ Oi as t glacial maximum — & Opre-industrial) : 

CAM3iso (e.g., Noone and Sturm, 2010; Pausata et al., 
2011a), ECHAM5-wiso (e.g., Werner et al., 2011), GISSE2- 
R (e.g., Schmidt et al., 2014; LeGrande and Schmidt, 2008, 
2009), IsoGSM (e.g., Yoshimura et al., 2003) and LMDZ4 
(e.g., Risi et al., 2010a). ECHAM5-wiso and IsoGSM 
outputs are for modern climate rather than pre-industrial 
conditions; however, the difference between the isotopic 
composition of pre-industrial and modern climate are 
expectedly small compared to late-glacial to late-Holocene 
<5 18 0 shifts. An offset factor was applied to simulated mean 
seawater <$ 18 0 in all five models (Table SI in the Supple- 
ment) to account for known glacial-interglacial changes 
to seawater <$ 18 0 (Emiliani, 1955; Dansgaard and Tauber, 
1969; Schrag et al., 1996, 2002). Possible spatial differences 
in seawater <S 18 0 changes from the last glacial maximum to 
the pre-industrial time period are not incorporated into simu- 
lations with prescribed sea surface temperatures (CAM3iso, 
ECHAM5-wiso, IsoGSM, LMDZ4) but are simulated by 
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the coupled ocean-atmosphere simulation of GISSE2-R 
(Table SI). GISSE2-R was submitted to the CMIP5 archive 
and participated in PMIP3. LMDZ4 was submitted to the 
CMIP3 archive. ECHAM5 and CAM3iso did not participate 
in CMIP5, while IsoGSM uses different boundary conditions 
than proposed for CMIP5 (Yoshimura et al., 2008). The 
five models span a range of spatio-temporal resolutions 
and isotopic/atmospheric parameterizations described in 
detail in the above references. A selection of the inter-model 
similarities and differences are summarized in Table S 1 . 

For clarity, empirical A 18 Oi at e-gi ac ial values that are 
based on measured isotope contents of groundwater, 
speleothem, ground ice or ice core records are referred 
to herein as measured A 18 Oi at e-gi ac iai; simulated pre- 
cipitation isotope compositions obtained from general 
circulation model results are referred to as simulated 
A 18 Oi a te-giacial- We acknowledge that the general circula- 
tion models explicitly analyse the last glacial maximum 
and the pre-industrial climate conditions (i.e., simulated 

A Oi a te-glacial = & Oi as t glacial maximum — ^ Op re -industrial)> 

whereas proxy record reconstructions of A 18 Oi at e-giacial 
integrate hydroclimatology over multi-millennial timescales 
that are different from the model simulations. 

3 Results and discussion 

3.1 Measured A 18 0| a te-giaciai values 

Measured groundwater (n = 59), speleothem (n = 15) 
and ice core (n = 12) A 18 Oi a t e - g iaciai values are pre- 
sented in Fig. 2 (references presented in the Supplement). 
Measured A 18 Oi at e-giaciai values range from —l.l%c 
(i.e., 5^ 8 0i a te-glacial < ^^Oi a te-Holocene) to +1.7 %c (i.e., 
5 18 Oiate-giaciai>5 18 Oiate-Hoioeene)- Three-quarters of the 
compiled records have negative measured A 18 Oi a t e -giaciai 
values and one-quarter of compiled records have positive 
measured A 18 Oi at e-gi ac iai values. Most groundwater-based 
late-glacial to late-Holocene shifts fall along <$ 2 H/<$ 18 0 
slopes of ~ 8 (Fig. S58 in the Supplement), suggesting that 
most groundwaters record temporal shifts to precipitation 
isotope contents rather than to soil evaporation isotope 
effects (see Evaristo et al., 2015). More than 80 % of records 
with positive measured A 18 Oi at e-gi ac iai values are located 
within 35° of the equator and within 400 km of the nearest 
coastline (e.g., Bangladesh A 18 Oi ate -giaciai of +1.5 %c, less 
than 300 km from the coast; Figs. 2-4). In comparison, 
negative measured A 18 Oi ate -giacial values are found in both 
coastal regions and farther inland. Negative measured 
A 18 Oi a te-giaciai values of the greatest magnitude are located 
at high latitudes (e.g., northwestern Canada, latitude 64° N: 
A 1 8 Oi a te-giaciai of — 5.5%c; northern Russia latitude 72° N: 
— 5.4 %c) and far from coastlines (e.g., Hungary: — 3.1 %c, 
~ 500 km from Atlantic Ocean; Peru: — 6.3 %c, ~ 2000 km 
from Atlantic Ocean, the modem moisture source to Peru; 
Garreaud et al., 2009). Greenland and Antarctic ice cores 
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Figure 2. Meteoric water <$ 18 0 change from the late-glacial 
(20000 to ~ 50 000 years ago) to the late-Holocene (within 
past ~5000 years; average A 18 Oi ate _ g i ac i a i values shown, 
where A^ 8 Oi a t e _gi ac i a i = <5^0i a t e _gi ac i a i — £^Oi a t e _H 0 i ocene ). The 
low temporal resolution of groundwater records means that <$ 18 0 
variations within each time period are smoothed and likely repre- 
sent unequal temporal weighting. References for measured mete- 
oric water <$ 18 0 changes for ice cores, groundwater and cave calcite 
are presented in the Supplement. 


have negative measured A 18 Oi a t e - g iaciai values that are of 
greater magnitude than non-polar measured A 18 Oi a t e -glacial 
values (Antarctic and Greenland A 18 Oi ate -giaciai values range 
from —3.6 to — 7.1 %o; Fig. 3). 

Our synthesis shows that measured A 18 Oi a t e -giaciai val- 
ues in the tropics are closer to 0 %c (i.e., no change) 
than A 18 0iate-gladal values at high latitudes and continen- 
tal interiors that generally have high magnitude, negative 
A 1 8 Oiate-giacial values. High magnitude, negative measured 
A 1 8 Oi a te-giaciai values are most common where present-day 
precipitation <S 18 0 values are at a minimum (e.g., Bowen 
and Wilkinson, 2002). This broad spatial pattern is consis- 
tent with the non-linear isotopic distillation of air masses un- 
dergoing progressive rainout (i.e., Rayleigh distillation). Be- 
cause seawater <5 18 0 values were ~ 1 %c higher- than-modem 
during the last glacial maximum (Schrag et al., 1996, 2002), 
our finding that the majority of measured A 18 Oi ate -giacial 
values are negative suggests that isotopic distillation of air 
masses was greater during the late-glacial than under present 
climate. This finding is consistent with land surface tem- 
perature reconstructions that show larger glacial-to-modem 
changes to land temperatures at high latitude and continental 
settings (Fig. 1; Annan and Hargreaves, 2013). Tropical ver- 
sus extratropical patterns of late-glacial/late-Holocene tem- 
perature change (Fig. la) are broadly similar to measured 
A 18 Oi a te-giaciai values (Fig. 3), where both temperature and 
isotope shifts are greater at high latitudes relative to the equa- 
tor. Therefore, it is possible that the larger late-glacial to late- 
Holocene temperature shifts at the poles relative to the equa- 
tor may have served to amplify the non-linear, Rayleigh rela- 
tionship describing the heavy isotope depletion of air masses 
undergoing progressive rainout during transport from lower 
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O Groundwater O Ice core A Cave 



CAM3iso ,-'ECHAM5 - 'GISSE2-R 
- - IsoGSM ^ LMDZ4 


Figure 3. Latitudinal variations of A 18 Oi ate _ g i ac i a i values of 
groundwater (circles, each circle is one aquifer), ice cores (dia- 
monds) and cave calcite (i.e., triangles; where A 18 Oi ate _ g i ac i a i = 
£ 1 8 Oiate-giaciai — 5l8 °late-Holocene)- Dashed lines mark 10° zonal 
mean simulated A 18 Oi ate _ g i ac i a i values from five different gen- 
eral circulation models: CAM3iso, ECHAM5-wiso, GISSE2-R, 
IsoGSM and LMDZ4 (Yoshimura et al., 2003; Legrande and 
Schmidt, 2008, 2009; Risi et al., 2010a; Noone and Sturm, 2010; 
Pausata et al., 201 la; Werner et al., 201 1). 


to higher latitudes. Further, the late-glacial was characterized 
by (i) lower-than-modern atmospheric temperatures with 
larger coastal-inland gradients, and (ii) lower-than-modern 
eustatic sea level leading to longer overland atmospheric 
transport distances. Each of these late-glacial/late-Holocene 
changes favours stronger-than-modern isotopic distillation 
of air masses transported inland from the coast during the 
late-glacial (Dansgaard, 1964; Rozanski, 1993; Winnick et 
al., 2014), potentially contributing to the broad, global ob- 
servation that most (77 %) 5 18 Oi ate -Hoiocene values exceed 
<5 18 0iate-glacial values on continents. 

Pairings of groundwater and speleothem records are 
available within ~ 500 km of one another in the south- 
western USA, central China and Israel. Southwestern 
USA speleothem and groundwater records ~ 400 km apart 
show similar A 18 Oi ate -giacial values, with San Juan Basin 
groundwaters having a measured A 18 Oi ate -giacial value of 
— 2.5±1.0%c (Phillips et al., 1986) and speleothems 
~ 400 km to the south having measured A 18 Oi a t e -giacial 
values of —3.0 ±1.2 and —3.4 ±0.4 (Asmerom et al., 
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2010; Wagner et al., 2010). Central China speleothem and 
groundwater records ~ 200 km apart overlap within uncer- 
tainty margins (i.e., A 18 Oi a t e -glacial values of —1.1 zb 1.7 
and +0.3 + 2.1 %c>; Cai et al., 2010). Israeli speleothem 
and groundwater records ~ 100 km apart have different 
measured A 18 Oi a t e -giadai values. Two Israeli groundwater 
A 1 8 Oi ate -gi a ciai records were compiled; the coastal Israeli 
aquifer has a A 18 Oi ate -gi ac ial value of +0.3 ± 0.4 %c (Yechieli 
et al., 2009), whereas groundwater of the Dead Sea Rift Val- 
ley has a A 18 Oi ate -gi a ciai value of —1.8 ±0 .6%c (Burg et al., 
2013). Speleothem records have A 18 Oi ate -gi a ciai values close 
to +1 %c (Frumkin et al., 1999; Bar-Matthews et al., 2003). 
In northern Turkey, speleothem and groundwater separated 
by ~ 150 km have measured A 18 Oi at e- g iacial values that dif- 
fer by ~3%c (speleothem A 18 Oi ate -giaciai — 5.7 ±0A%c ver- 
sus groundwater A 18 Oi ate -giacial of —2.8 ± 1.0 %c\ Fleitmann 
et al., 2009; Arslan et al., 2013, 2015). While the locations of 
the groundwater and speleothem records differ, the compiled 
data suggest that groundwater and speleothem A 18 Oi ate -giacial 
values may capture different A 18 Oi ate -gi a cial values under 
similar climate conditions. 

A number of potential processes could bias the preser- 
vation of precipitation isotope composition in ice core, 
speleothem or groundwater archives (Wang et al., 2001; 
Thompson et al., 2006; Edmunds, 2009). For example, 
groundwater and speleothem archives preserve only the iso- 
tope record of precipitation that traverses the vadose zone. 
Recent global analyses of paired precipitation-groundwater 
isotope compositions show that winter (extratropics) and wet 
season (tropics) precipitation contributes disproportionately 
to recharge (Jasechko et al., 2014), meaning that paleocli- 
mate records may be more sensitive to changes to winter and 
wet seasons than summer or dry season (Vogel et al., 1963; 
Simpson et al., 1972; Grabczak et al., 1984; Harrington et al., 
2002; Jones et al., 2000; Darling, 2004; Partin et al., 2012). 
Similarly, groundwater isotope records are unlikely to repre- 
sent constant and continuous recharge fluxes during the late- 
Holocene or the late-glacial (McIntosh et al., 2012). Mod- 
em groundwater recharge fluxes are highest in humid cli- 
mates (Wada et al., 2010). Groundwater <$ 18 0 records only 
represent precipitation that recharges aquifers, meaning that 
groundwater-based A 18 Oi ate -gi a ciai values could be biased to 
subintervals (e.g., interstadials, pluvial periods) within the 
late-Holocene and late-glacial intervals when recharge fluxes 
were at local maxima. Speleothem records may be further 
complicated by processes impacting the timing of calcite pre- 
cipitation. Recent modelling suggests that calcite precipita- 
tion in caves located outside of the tropics is greatest dur- 
ing the cool season and reduced during summer months 
due to changes in ventilation, meaning that higher latitude 
speleothems record oxygen isotope compositions biased to 
cool season climate change (James et al., 2015). Other recent 
work suggests that speleothem <5 ls O data may be impacted 
by disequilibrium isotope effects (Asrat et al., 2008; Daeron 
et al., 2011; Kluge and Affek, 2012; Kluge et al., 2013) 
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or by partial evaporation of drip waters resulting in 18 0- 
enrichment (e.g., Cuthbert et al., 2014a) and greater frac- 
tionation due to evaporative cooling (Cuthbert et al., 2014b), 
potentially obscuring the preservation of primary precipita- 
tion isotope contents in the speleothem record. Compiled ice 
core records may have been influenced by post-depositional 
exchanges of ice with atmospheric vapour (Steen-Larsen et 
al., 2014). The impact of atmospheric vapour exchanges on 
ice core isotope records remains poorly understood. Potential 
biases in the preservation of precipitation <$ 18 0 differ among 
groundwater, glacial ice, and speleothem records, meaning 
that co-located records of differing record-type may pre- 
serve different A 18 Oi ate -gi a cial values under similar climate 
conditions. Finally, all proxy records may be impacted by 
past changes in the seasonality of precipitation, which can 
substantially impact annual precipitation <$ 18 0 values (e.g., 
Werner et al., 2000). 

We cannot rule out the possibility that changes in sea- 
sonal biases of proxy record preservation occurred between 
the late-glacial and the late-Holocene and have impacted 
measured A 18 Oi ate -gi a cial values. Further, the chronologies 
of groundwaters and ice core records have uncertainties 
on the order of thousands of years, meaning that the time 
intervals used to calculate measured A 18 Oi ate -gi a ciai values 
may be inaccurate. However, the plateauing of isotope con- 
tent observed in most regional aquifers for 0-5000 years 
before present and for >20000 years before present sup- 
ports our interpreting these data as records of late-glacial 
to late-Holocene isotopic shifts (see figures in the Supple- 
ment). Notwithstanding potential <$ 18 0 preservation biases 
and chronology uncertainties, the global data synthesized 
here show patterns consistent with the enhanced distillation 
of advected air masses originating as (sub)tropical ocean 
evaporate and undergoing progressive, poleward rainout un- 
der cooler-than-modern late-glacial temperatures. 

3.2 Simulated A^ ^O^e-giaciai values 

Simulated precipitation A 18 Oi ate -gi a cial values from five gen- 
eral circulation models are presented in Fig. 5. At least 
four of the five models agree on the sign of simulated 
A 1 8 Oi a te-gi ac iai values - that is values consistently above or 
consistently below zero - for 68.8 % of grid cells cover- 
ing Earth’s surface (68.7 % of over-ocean areas and 68.9 % 
of over land areas; multi-model calculation completed using 
three of four models as a threshold at high-latitudes where 
IsoGSM data were unavailable). Simulated A 18 Oi ate -giaciai 
values are consistently negative over the North Atlantic 
Ocean and the Fennoscandian and Laurentide ice sheets 
and consistently positive over most of the tropical oceans, 
whereas poorer agreement is found over tropical land sur- 
faces. The negative simulated A 18 Oi ate -gi ac ial values over the 
Northern Hemisphere ice sheets and North Atlantic are likely 
driven by the difference in ice sheet topography and sea 
ice cover, between the late-glacial and pre-industrial cli- 
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O Groundwater A Cave 



Figure 4. Measured A 18 Oi ate _ g i ac i a i value variability with dis- 
tance to the nearest coast (A 18 Oi ate _ g i acia i = <$ 18 Oi ate _ g i acia i - 
<5 1 8 Oiate-Holocene) • Tropical and subtropical locations are shown in 
deep blue (<35° absolute latitude), extra-tropical sites are shown 
in light grey (> 35° absolute latitude). The shape of each point cor- 
responds to groundwater and ground ice (circles) or cave calcite 
(i.e., speleothems; triangles). Error bars mark one standard devia- 
tion from the mean. 


mate. The late-glacial to late-Holocene change in ice sheet 
topography and sea ice cover impacted surface tempera- 
tures, which were more than ~20°C cooler over most of 
present-day Canada during the last glacial maximum (Fig. 1). 
Cooler temperatures in conjunction with ice sheet topog- 
raphy (> 3000 m elevations; e.g., Peltier, 1994) enhanced 
Rayleigh distillation for air masses transecting Northern 
Hemisphere ice sheets, as evidenced by systematically low 
measured and simulated <$ 18 Oi ate -giacial values in these re- 
gions (Figs. 2, 3 and 5). 

A comparison of simulated A 18 0i a te-gladal values over 
tropical Africa, South America and Oceania shows inter- 
model disagreement (Fig. 5). Different tropical simulated 
A 1 8 Oiate-giadal values among the models reflect the different 
isotopic parameterizations, inter-model spread in simulated 
precipitation rates, and seawater <$ 18 0 specifications used in 
each model (Supplement). Inter-model spread in simulated 
A 1 8 Oi a te-giacial values in some regions highlights the impor- 
tance of this global synthesis of measured A 18 Oi ate -giacial 
values as a constraint for isotope-enabled climate simula- 
tions. Another potential source for the model disagreement 
is introduced by the different ice- sheet topography used in 
each model. CAM3Iso, IsoGSM and LMDZ4 used Ice 5G 
(Peltier, 1994) as advised for PMIP2 (Braconnot et al., 2007), 
whereas the GISSE2 replaces Ice 5G Laurentide ice with that 
of Licciardi et al. (1999) and ECHAM5-wiso uses ice topog- 
raphy from PMIP3 (Braconnot et al., 2007, 2012; PMIP3 fol- 
lows ice sheet topography blended from multiple ice sheet re- 
constructions: Argus and Peltier, 2010; Toscano et al., 201 1). 
Ice sheet topography is an important driver of simulated tem- 
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perature, precipitation and atmospheric circulation during the 
last glacial maximum (e.g., Justino et al., 2005; Pausata et 
al., 2011b; Ullman et al., 2014). Therefore, it is likely that 
inter-model differences in paleo-ice sheet topographies im- 
pacts atmospheric circulation and thus high latitude simu- 
lated A 1 8 Oi a te-giacial values reported in this study (Fig. 5). 

Differences in the specification of initial seawater <$ 18 0 
may also lead to inter-model differences in simulated 
A 1 8 Oi a te-giaciai values. Seawater <$ 18 0 is set to be globally 
homogenous in CAM3Iso, IsoGSM and LMDZ4, and het- 
erogeneous in ECHAM5-wiso (using modern gridded sea- 
water <$ 18 0 heterogeneity of LeGrande and Schmidt, 2006) 
and GISSE2-R (coupled atmosphere-ocean model; seawa- 
ter <S 18 0 is calculated by the ocean model). Including sur- 
face ocean <$ 18 0 heterogeneities in model simulations im- 
pacts land precipitation <$ 18 0 by up to ~ 1.5 %c relative to 
simulations with homogenous seawater <$ 18 0 (LeGrande and 
Schmidt, 2006). However, different seawater <$ 18 0 specifica- 
tions cannot account for all inter-model differences in simu- 
lated A 1 8 Oi a te-glacial values. 

The models also show deficiencies in simulating measured 
A 18 Oi a te-giacial values in the tropics, particularly over trop- 
ical Africa. This finding could, in part, relate to the high 
sensitivity of precipitation 8 18 O to convective parameteriza- 
tions (Lee et al., 2009; Field et al., 2014), although future 
research is required to test this. Another reason may be that 
the measured A 18 Oi a t e -giacial integrates the hydroclimatologi- 
cal signal over multi-millennial timescales, whereas the sim- 
ulated A 1 8 Oi a te-giacial values explicitly explore last glacial 
maximum and pre-industrial/present-day climate conditions. 
The smeared temporal resolution of groundwater-based mea- 
sured A 1 8 Oi a te-giacial values due to storage and mixing in the 
aquifer precludes an ideal comparison of measured versus 
simulated A 18 Oi ate -giacial values. Further, as previously dis- 
cussed in Sect. 3.1, the measured A 18 Oi a t e -giaciai values are 
susceptible to a number of potential biases that may obscure 
the magnitude and direction of late-glacial to late-Holocene 
precipitation <$ 18 0 changes. Notwithstanding, models cor- 
rectly simulate the sign of measured A 18 Oi ate -giacial values 
(i.e., positive or negative) in the extratropics more frequently 
than in the tropics. Better agreement in the sign of sim- 
ulated versus measured A 18 Oi ate -giaciai values in the extra- 
tropics compared to the tropics is likely linked to the substan- 
tial changes to extra-tropical ice-sheet topography and sea- 
ice cover between the two climate states in northern North 
America and Europe. Substantial changes to Northern Hemi- 
sphere ice volumes between the late-glacial and the late- 
Holocene likely enhanced upwind distillation of air masses 
leading to high-magnitude, negative A 18 Oi ate -giacial values 
that are well captured by the climate simulations. However, 
simulated A 18 Oi ate -giaciai values over Antarctica and Green- 
land show large inter-model spread, suggesting that model- 
based interpretations of polar ice core records may vary 
widely among different atmospheric models. 
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Figure 5. Simulated precipitation <5 18 0 differences between the last glacial maximum and pre-industrial time periods (i.e., 
<5 18 0last glacial maximum - <5 18 0pre-industrial) f rom ^ ve general circulation models: CAM3iso, ECHAM5-wiso, GISSE2-R, IsoGSM and 
LMDZ4 (Yoshimura et al., 2003; Legrande and Schmidt, 2008, 2009; Risi et al., 2010a; Noone and Sturm, 2010; Pausata et al., 2011a; 
Werner et al., 2011). Circles (groundwater), triangles (speleothems) and diamonds (ice cores) show measured A 18 Oi ate _ g i ac i a i values from 
paleoclimate proxy records (Fig. 1, original data presented in Tables S2-S5). The panel entitled “Composite” shows the multi-model ensem- 
ble median simulated A 18 Oi ate _ g i ac i a i value where at least four of the five models agree on the sign of simulated A 18 Oi ate _ g i ac i a i values (i.e., 
positive or negative; all five model simulations of 5 18 0] ast g i ac i a i maximum - ^ 18 Opre-industrial were used to calculate multi-model median 
shown in “Composite”). 


3.3 Regional measured and simulated A 18 0| a te-giaciai 
values 


3.3.1 Australia and Oceania 

Measured A 18 Oi a t e -giadai values from Australia and Oceania 
fall between —1 and +1 %c (Fig. 2). Australian climate dur- 
ing the last glacial time period was more arid (Nanson et 
al., 1992), dustier (Chen et al., 1993) and cooler (Miller et 
al., 1997) than present day. Simulated A 18 Oi ate -gi a ciai values 
across Australia are variable among the five models. Mea- 
sured A 1 8 Oi ate -gi a ciai values across Oceania have been at- 
tributed to temporal changes in the strength of monsoons and 
convective rains (Aggarwal et al., 2004; Partin et al., 2007; 
Williams et al., 2010) potentially impacted by late-glacial to 
late-Holocene shifts in the position of the intertropical con- 
vergence zone (Lewis et al., 2010, 2011). 


3.3.2 Southeast Asia 

Measured A 18 Oi ate -gi a ciai values from southeastern Asia 
range from —2.3 to +1.7 %c . The highest regional mea- 
sured A 1 8 Oi a te-giacial values are found in Bangladesh (mea- 
sured A 1 8 Oi a te-giaciai of +1.5 + 1.3 %o\ Aggarwal et al., 
2000) and in central and southeastern China (measured 
A 1 8 Oi a te-giaciai of +0.3 to +1.7 %o\ Wang et al., 2001; Yuan 
et al., 2004; Dykoski et al., 2005; Cai et al., 2010; Yang et 
al., 2010). General circulation models have positive simu- 
lated A 1 8 Oi a te-giaciai values near to the Chinese coasts, but 
are more variable across western and northern China (Fig. 5). 
Chinese speleothem records show near-zero or positive mea- 
sured A 1 8 Oi ate -gi a cial values interpreted to reflect the reduced 
strength of the East Asian (Wang et al., 2001; Dykoski 
et al., 2005; Cosford et al., 2008) or Indian monsoons 
(Pausata et al., 2011a). Further research suggests that Chi- 
nese speleothem 5 18 0 variations reflect changes to regional 
moisture sources and the intensity or provenance of atmo- 
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spheric transport pathways (LeGrande and Schmidt, 2009; 
Dayem et al., 2010; Lewis et al., 2010; Maher and Thomp- 
son, 2012; Caley et al., 2014b; Tan, 2014). 

North China Plain groundwaters have high-magnitude, 
negative A 18 Oi a t e -giaciai values (measured A 18 Oi a t e -giaciai of 
—2.3 zb 0.6 %o\ Chen et al., 2003) compared to coastal, more 
southerly counterparts. Combining the negative measured 
A 18 Oi ate -giacial in northern China (Chen et al., 2003; Ma et al., 
2008; Currell et al., 2012; Li et al., 2015) with the positive 
measured A 18 Oi ate -gi a cial values in central and southeastern 
China (Wang et al., 2001; Yuan et al., 2004; Dykoski et al., 
2005; Cai et al., 2010; Yang et al., 2010) reveals a south-to- 
north decrease from positive (south) to negative (north) mea- 
sured A 1 8 Oi a te-giacial values (Figs. 2 and 6). Previous stud- 
ies of modern precipitation have identified increasing pre- 
cipitation <$ 18 0 values from the coast to inland China during 
the wet season, sharply contrasting spatial patterns expected 
from Rayleigh distillation (Araguas-Araguas et al., 1998). A 
more recent work suggests that low wet- season precipitation 
<5 18 0 values over southern China are controlled by the de- 
flection of westerlies around the Tibetan Plateau, whereas 
precipitation <$ 18 0 values over northern China are controlled 
by local-scale rainfall and below-cloud raindrop evapora- 
tion (Lee et al., 2012). Therefore, measured A 18 Oi ate -giaciai 
values from southern China may reflect changes to atmo- 
spheric circulation at broader spatial scales, whereas mea- 
sured A 18 Oiate-giacial values from northern China may indi- 
cate changes to more localized atmospheric conditions im- 
pacting processes such as raindrop evaporation in addition to 
meso- and synoptic-scale circulation changes. 

3.3.3 Africa 

Measured A 18 Oi ate -gi a cial values from Africa range from 
—2.9 to +0.1 %o (Figs. 2 and 6). Sixteen of 17 measured 
A 18 Oi a te-giacial values from Africa are negative. Near-zero 
measured A 18 Oi ate -giaciai values are generally found near to 
coasts (e.g., Senegal A 18 Oi a t e - g iaciai of +0.1 ±0.8 %c\ Ma- 
dioune et al., 2014), whereas higher magnitude, negative 
measured A 18 Oi ate -giacial values in Africa are found farther 
inland (e.g., Niger A 18 Oi a t e - g iacial values of —2.3 ±2.0 and 
— 2.9±0.9%c: ~ 800 km from the Atlantic coast). General 
circulation model A 18 Oi ate -giaciai values show poor agree- 
ment with measured A 18 Oi a te-glacial over tropical Africa com- 
pared to model-measured comparisons for Europe and North 
America (Fig. 5), with positive simulated A 18 Oi a te-glacial 
values predicted over large parts of Africa where negative 
A 18 Oi a te-giaciai values are measured. Figure 5 shows that 
Africa has the largest inter-model and model-measurement 
disagreements in the sign of A 18 Oi ate -giacial values of the con- 
tinents. 

Northern African hydrological processes are influenced by 
interlinked controls such as meridional shifts in the posi- 
tion of the intertropical convergence zone (Arbuszewski et 
al., 2013) and the strength of Atlantic meridional overturning 



-5.5 -4.5 -3.5 -2.5 -1.5 -0.5 +0.5 +1.5 
Late-glacial 5 18 0 minus late-Holocene 5 18 0 (%o SMOW) 


Figure 6. Regional proxy record A 18 Oi ate _ g i ac i a i values for 
(a) southeastern Asia, (b) Africa, (c) Europe, and (d) the 
contiguous United States of America (where A 18 Oi ate _ g i ac i a i = 
5 1 8 Oi a te-giacial — ^Olate-Holocene)- The multi-model ensemble 
median simulated A 18 Oi ate _gi ac jai value is shown as a grid (0.5 de- 
gree smoothing). Groundwater records are represented by circles, 
speleothems by triangles, and ice cores by diamonds, labels show 
measured A 18 Oi ate _ g i ac i a i values for each individual record. 


circulation (Mulitza et al., 2008). Paleowater chemistry indi- 
cates that northern Africa was at least 2 °C cooler than today 
(Guendouz et al., 1998) and that westerly moisture transport 
was stronger than the present during the late-glacial (Sultan 
et al., 1997; Abouelmagd et al., 2012). 

Tropical Africa was 2 to 4°C cooler and more arid than 
present day at the last glacial maximum (Powers et al., 
2005; Tierney et al., 2008). Early- and late-Holocene rain- 
fall and isotope compositions were highly variable across 
Africa (Tierney et al., 2008, 2013; SchefuB et al., 2011; Otto- 
Bliesner et al., 2014). Tropical African rainfall originates 
from both Indian and Atlantic sources, with Atlantic- sourced 
moisture travelling across the Congo rainforest (Levin et al., 
2009). Lower- than-modern continental moisture recycling 
during the late-glacial may partially explain negative mea- 
sured A 1 8 Oi ate -gi a cial values across some regions of inland 
tropical Africa (e.g., Risi et al., 2013). Negative measured 
A 1 8 Oi a te-gi ac ial values in tropical Africa could also be in- 
terpreted to reflect higher-than-modern upwind rainout dur- 
ing the late-glacial (see Risi et al., 2008, 2010b; Lee et al., 
2009; Scholl et al., 2009; Lekshmy et al., 2014; Samuels- 
Crow et al., 2014); however, this explanation necessitates 
stronger-than-modern convection during the late-glacial, an 
explanation that would contradict the established cooler- 
than-modern land surface temperatures. Therefore, changes 
to atmospheric transport distances and vapour origins are 
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more likely responsible for negative measured A 18 Oi a t e -glacial 
values across tropical Africa (Lewis et al., 2010). 


3.3.4 Europe and the Mediterranean 

Measured A 18 Oi at e-giaciai values across Europe, the Middle- 
East and the eastern Mediterranean range from —5.7 to 
+ 1.3 %o. Eighty percent of measured A 18 Oi ate -gi a cial values 
across these regions are negative. All five general circula- 
tion models agree on negative simulated A 18 Oi ate -gi a ciai val- 
ues across Europe, consistent with the negative measured 
A 1 8 Oi ate -gi a cial values across the majority of Europe. Mea- 
sured A 1 8 Oi a te-giacial values are generally higher in west- 
ern Europe (0.0 to — 1.0 %c in Portugal, the United King- 
dom and France) than in eastern Europe (—1.0 to —5.1%c 
in Poland, Hungary and Turkey; Stute and Deak, 1989; Le 
Gal La Salle et al., 1996; Darling et al., 1997; Barbecot et 
al., 2000; Zuber et al., 2004; Galego Fernandes and Carreira, 
2008; Celle-Jeanton et al., 2009; Varsanyi et al., 2011; Sam- 
borska et al., 2013; Arslan et al., 2013). This spatial pattern 
of A 1 8 Oi a te- g i ac ial values is consistent with enhanced isotopic 
distillation of westerlies during the late-glacial due to cooler- 
than-modern final condensation temperatures. 

High magnitude, negative measured A 18 Oi ate -gi a cial values 
are located in Turkey and Georgia south and east of the Black 
Sea (—2.8 ±1.0 to — 5.7±0.4%c; Fleitmann et al., 2009; 
Arslan et al., 2013; Melikadze et al., 2014). Westerly air 
mass trajectories distal to the Fennoscandian ice sheet to- 
pography may not have changed considerably since the late- 
glacial over western and central Europe (Rozanski, 1985; 
Loosli et al., 2001). Therefore, higher, near-zero measured 
A 1 8 Oi a te-gi ac iai values in western Europe and lower, negative 
measured A 18 Oi ate -gi a cial values in eastern Europe indicate 
enhanced distillation of advected air masses during the late- 
glacial relative to the late-Holocene. 

Changes to freeze-thaw conditions of the ground sur- 
face between the latter half of the last glacial time period 
and the modern climates may have impacted the seasonal- 
ity of the fraction of precipitation recharging aquifers and 
thus A 1 8 Oi ate - g i a ciai (Darling, 2004, 2011; Jasechko et al., 
2014). Geomorphic evidence suggests permafrost covered 
portions of Hungary at the last glacial maximum, suggest- 
ing that land temperatures may have been up to 15 °C cooler 
than present day (Fabian et al., 2014), a larger late-glacial 
to late-Holocene temperature shift than earlier, noble gas- 
based reconstructions (5-7 °C; Deak et al., 1987). European 
pollen records indicate that northern Europe was tundra- 
like and that southern Europe was semi-arid during the last 
glacial maximum (Harrison and Prentice, 2003; Clark et al., 
2012). The European late-glacial to late-Holocene transition 
from semi-arid deserts to temperate forests could have low- 
ered A 1 8 Oi a te-gi ac iai values as groundwater recharge ratios 
transitioned from more extreme winter-biased (e.g., semi- 
arid lands during the late-glacial) to less extreme winter- 
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biased groundwater recharge ratios (e.g., forests during late- 
Holocene; Jasechko et al., 2014). 

3.3.5 South America 

Measured A 18 Oi ate -gi a ciai values across South America range 
from —6.3 to ±0.6 %c (Figs. 2 and 6). The highest- 
magnitude, negative measured A 18 Oi ate -gi a ciai values are 
found in Andean ice cores (A 18 Oi ate -giacial of —4.6 ±1.0 and 
—6.3 ± 1.3; Thompson et al., 1995, 1998). Here the impor- 
tance of upstream convection upon modem Andean precipi- 
tation 5 18 0 has been highlighted at inter- annual (Hoffmann 
et al., 2003; Vuille and Werner, 2005), seasonal (Vimeux et 
al., 2005; Samuels-Crow et al., 2014) and daily timescales 
(Vimeux et al., 2011). It is therefore possible that upstream 
convection controls past changes to Andean precipitation iso- 
tope compositions recorded in ice cores. 

The measured groundwater A 18 Oi ate -gi a cial value located in 
eastern Brazil is —2.7 ± 1.3 %o (Salati et al., 1974). Eastern 
Brazil was 5 °C cooler than today during the latter half of 
the last glacial period (Stute et al., 1995b). Four of the five 
general circulation models simulate positive A 18 Oi at e-gi ac ial 
values across eastern Brazil (Fig. 5), highlighting a differ- 
ence between simulated and measured A 18 Oi ate -giaciai values 
in parts of the tropics. The negative measured A 18 Oi ate -giacial 
value in eastern Brazil has been previously interpreted to re- 
flect higher-than-modem precipitation during the last glacial 
time period (Salati et al., 1974). Lewis et al. (2010) show 
that localized rainfall governs precipitation <5 18 0 in eastern 
Brazil. Modern precipitation 5 18 0 values are lowest in east- 
ern Brazil when precipitation rates are at a maximum. Ex- 
tending Lewis et al.’s interpretation linking local precipita- 
tion amount to precipitation <$ 18 0 would suggest that the neg- 
ative measured A 18 Oi ate -gi a cial value found in eastern Brazil 
may indeed record wetter-than-modem conditions during the 
late-glacial as proposed by Salati et al. (1974). Further, dis- 
agreement between measured and simulated A 18 Oi at e-gi ac iai 
in eastern Brazil highlights the need to critically evaluate cli- 
mate model performance in regions where the precipitation 
amount is closely correlated with precipitation <$ 18 0. 

3.3.6 North America 

Measured A 18 Oi ate -giacial from North American proxy 
records range from —5.5 to ±1.0%c. Canadian records of 
groundwater recharge that took place beneath the Lauren- 
tide ice sheet are not included in this synthesis (“subglacial 
recharge”; Grasby and Chen, 2005; Ferguson et al., 2007; 
McIntosh et al., 2012; Ferguson and Jasechko, 2015). These 
records were excluded because the subglacial meltwaters 
that recharged aquifers likely reflect precipitation that fell 
elsewhere on the paleo-ice sheet, potentially complicating 
the comparison of groundwater isotope compositions for the 
late-Holocene and last glacial time period. 
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Measured A 18 Oi a t e -giadai values along the USA east coast 
show the highest, positive values in Georgia (latitude: 32° N; 
measured A 18 Oi a t e -giaciai of +1.0 %c\ Clark et al., 1997), 
decreasing northward to near-zero measured A 18 Oi ate -gi a cial 
values in coastal Maryland (latitude 39° N; measured 
A 1 8 Oi ate -gi a cial of — 0.1 + 0.4 %c; Aeschbach-Hertig et al., 
2002). Decreasing A 18 Oi ate -gi a cial values with increasing lat- 
itude along the USA east coast may be explained in part 
by the isotopic distillation of air masses advected north- 
ward from the subtropics under cooler-than-modem final at- 
mospheric condensation temperatures. Indeed, paleoclimate 
records indicate that Maryland was more arid and as much 
as 9-12 °C cooler during the late-glacial relative to the late- 
Holocene (Purdy et al., 1996; Aeschbach-Hertig et al., 2002; 
Plummer et al., 2012). In addition to temperature change, 
late-glacial precipitation isotope compositions along east- 
ern USA coastline were likely impacted by the lower-than- 
modern late-glacial sea levels, which changed overland at- 
mospheric transport distances between the late-glacial and 
late-Holocene (Clark et al., 1997; Aeschbach-Hertig et al., 
2002; Tharammal et al., 2013). 

Measured A 18 Oi at e-gi ac ial values in the central and south- 
western USA have the highest magnitude, negative measured 
A 1 8 Oi ate -gi a cial values of temperate North America, ranging 
from —1.0 to — 3A%c . Central and southwestern USA mea- 
sured A 1 8 Oi a te-giaciai values contrast the positive measured 
A 18 Oi a te-giadai values found along the eastern USA coast at 
similar latitudes. Consistently negative A 18 Oi at e-giacial values 
in central and southwest USA suggest that advected moisture 
to the region underwent greater upstream air mass distillation 
during the late-glacial than under modem climate. Pollen, 
vadose zone and groundwater records show that late-glacial 
southwestern USA was ~4°C cooler, had greater ground- 
water recharge fluxes, and had more widespread forests than 
present day (Stute et al., 1992, 1995a; Scanlon et al., 2003; 
Williams, 2003). Negative measured A 18 Oi at e-gi ac iai values 
found in the southwest USA have been ascribed to lower- 
than-modern summer precipitation (New Mexico, Phillips et 
al., 1986), latitudinal shifts in the positions of the polar jet 
stream and the intertropical convergence zone (New Mexico, 
Asmerom et al., 2010) and changes to over-ocean humid- 
ity, temperature or moisture sources (Idaho, Schlegel et al., 
2009). Wagner et al. (2010) interpret decreases to southwest- 
ern precipitation <$ 18 0 to reflect cooler and more-humid con- 
ditions. Extending this interpretation to negative measured 
A 1 8 Oi a te- g i ac ial values found across the southwestern USA 
values supports earlier conclusions that the region was cooler 
and more humid than today during the late-glacial, possi- 
bly linked to changes in air mass trajectories and moisture 
sources (Asmerom et al., 2010; Wagner et al., 2010). Sim- 
ulated A 1 8 Oi a te-gi ac iai values across North America closely 
match spatial patterns of measured A 18 Oi ate -gi a cial synthe- 
sized in this study. Strong, multi-model agreement with mea- 
sured A 1 8 Oi ate -gi a cial patterns supports continued application 
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of isotope-enabled general circulation models when inter- 
preting North American precipitation isotope proxy records. 

4 Conclusions 

While changes to the isotope content of precipitation be- 
tween the last glacial time period and more recent times 
has been widely documented, few studies have synthesized 
these dispersed data to explore the global patterns of <$ 18 0 
change driven by past shifts to regional climate. In this study 
we compile groundwater, speleothem, ice core and ground 
ice records of <$ 18 0 shifts between the late-glacial (20 to 
~ 50 thousand years ago) and the late-Holocene (within the 
past 5000 years). Late-glacial to late-Holocene <$ 18 0 shifts 
range from —7.1 %c (i.e., <^ 8 Oi a te-gl a cial <^ 8 Oi a te-Holocene) to 
+ 1.7 (i.e., S 18 Oiate-giaciai>5 18 Oiate-Hoiocene)- Aquifers with 
positive measured A 18 Oi ate -giaciai values (23% of records) 
are most common along the subtropical coasts. The majority 
(77 %) of measured A 18 Oi ate -giaciai values are negative, with 
the highest magnitude differences between 5 18 0i a t e - g iacial 
and 5 1 8 Oi a te-Holocene observed at high latitudes and far from 
coasts. This spatial pattern suggests that isotopic distillation 
of advected air masses was greater during the late-glacial 
than under present climate, likely due to the non-linear na- 
ture of Rayleigh distillation, accentuated by larger glacial- 
interglacial atmospheric temperature changes at the poles rel- 
ative to lower latitudes. Regionally divergent precipitation 
<$ 18 0 responses to the ~4°C of global warming occurring 
between the late-glacial and the late-Holocene suggest that 
continued monitoring of modern precipitation isotope con- 
tents may prove useful for detecting hydrologic changes due 
to ongoing, human-induced climate change. Future paleo- 
precipitation proxy record <$ 18 0 research can use these new 
global maps of A 18 Oi ate -gi ac ial records to target and prior- 
itize field sites. In the near term, a global compilation of 
large lake sediment isotope records that accounts for paleo- 
evaporative isotope effects could enhance spatial coverage of 
interglacial-glacial <5 18 0 shifts. 

General circulation models agree on the sign and magni- 
tude of terrestrial precipitation A 18 Oi ate -gi a cial values better 
in the extra-tropics than in the tropics. Differences in simu- 
lated precipitation isotope composition changes amongst the 
models might be linked to different parameterizations of sea- 
water <5 18 O, glacial topography and convective rainfall, how- 
ever, these hypotheses require further testing. Future model 
research should focus on quantifying the relative roles of 
inter-model spread in the simulated climate versus the iso- 
topic response to climate change on resulting simulated pre- 
cipitation <$ 18 0. This would provide guidelines to interpret 
model-data isotopic differences and to identify what aspects 
climate models have greatest difficulties capturing. 

The Supplement related to this article is available online 
at doi:10.5194/cp-ll-1375-2015-supplement. 
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